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ABSTRACT 

Models of the chemical evolution of our Galaxy are extended to include radial mi- 
gration of stars and flow of gas through the disc. The models track the production 
of both iron and a elements. A model is chosen that provides an excellent fit to the 
metallicity distribution of stars in the Geneva-Copenhagen survey (GCS) of the solar 
neighbourhood, and a good fit to the local Hess diagram. The model provides a good 
fit to the distribution of GCS stars in the age-metallicity plane although this plane 
was not used in the fitting process. Although this model's star- formation rate is mono- 
tonic declining, its disc naturally splits into an a-enhanced thick disc and a normal 
thin disc. In particular the model's distribution of stars in the ([0/Fe],[Fe/H]) plane 
resembles that of Galactic stars in displaying a ridge line for each disc. The thin-disc's 
ridge line is entirely due to stellar migration and there is the characteristic variation of 
stellar angular momentum along it that has been noted by Haywood in survey data. 
Radial mixing of stellar populations with high from inner regions of the disc to 
the solar neighbourhood provides a natural explanation of why measurements yield 
a steeper increase of az with age than predicted by theory. The metallicity gradient 
in the ISM is predicted to be steeper than in earlier models, but appears to be in 
good agreement with data for both our Galaxy and external galaxies. The models are 
inconsistent with a cutoff in the star-formation rate at low gas surface densities. The 
absolute magnitude of the disc is given as a function of time in several photometric 
bands, and radial colour profiles are plotted for representative times. 

Key virords: galaxies: abundances - galaxies: evolution - galaxies: ISM - galaxies: 
kinematics and dynamics - Galaxy: disc - solar neighbourhood 



1 INTRODUCTION 

Models of the chemical evolution of galaxies are key tools 
in the push to understand how galaxies formed and have 
evolved. Their application to our Galaxy is of particular im- 
portance both on account of the wealth of observational data 
that they can be required to reproduce, and on account of 
the inherent interest in deciphering the history of our envi- 
ronment. 

From the pione ering papers by Ivan den Berghl (Il962l ) 
and ISchmidtl (|l963 l) it has generally been assumed that a 
galaxy such as the Milky Way can be divided into concen- 
tric cylindrical an nuli, each of which evolves independently 
of the others (e.g. |Pagellll997l: IChiappini et all 1 19971 . l200ll : 



iNaab fc Ostrike j l2006l : IColavitti et all 120081 ). The contents 
of any given cylinder are initially gaseous and of extremely 
low or zero metallicity. Over time stars form in the cylinder 
and the more massive ones die, returning a mixture of heavy 
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elements to the remaining gas. The consequent increase in 
the metallicity of the gas and newly-formed stars is gener- 
ally moderated by an inflow of gas from intergalactic space, 
and, less often, by an outflow of supernova-heated gas. 

The cool, star-forming gas within any cylinder is as- 
sumed to be well mixed, so at any time it can be char- 
acterised by a metallicity Z{r,t), where r is the cylin- 
der's radius. Hence the stars formed within a given cylinder 
should have metallicities Z{r, tf) that are uniquely related to 
their time of formation, tf . Observations do not su bstanti- 
ate this prediction; in fact lEdvardsson et al.l (|l993l ) showed 
that solar-neighbourhood stars are widely distributed in the 
{t[, Z) plane - for a detailed discussion see iHavwood ( 2006h 
and Section [6.21 

The absence of an age-metallicity relation in the so- 
lar neighbourho od is naturally explained by radial mi- 
gration of stars (ISellwood fc Binnevl I2OO2I : iHavwoodl I2OO8I : 
iRoskar et al.ll2008bl ). It has been recognised for many years 
that scattering by spiral structure and molecular clouds 
gradually heats the stellar disc, moving stars onto ever more 
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eccentric and inclined orbits. Stars that are on eccentric or- 
bits clearly contribute to different cylindrical annuli at dif- 
ferent phases of their orbits, and thus tend to modify any 
radial gradient in the metallicities of newly formed stars. 
Moreover, scattering events also change the guiding centres 
of stellar orbits, so even a star on a circular orbit can be 
found at a differen t radiu s from that of its birth. In fact, 
ISellwood fc BinnevI (|2002l ) argued that the dominant effect 
of transient spiral structure is resonant scattering of stars 
across the structure's corotation resonance, so even a star 
that is still on a near-circul ar orbit may be far from its 
radius of birth. iRoskar et al. 1 iioOSa) showed that in a cos- 
mological simulation of galaxy formation that included both 
stars and gas, resonant scattering at corotation caused stars 
to move outwards and gas inwards, with the result that the 
stellar disc extended beyond the outer limit of star forma- 
tion; the outer disc was entirely populated by stars that 
had formed much further in and yet were still on nearly 
circular orbits. This simu lation confirmed the conjecture of 
ISellwood fc BinnevI l|2002l ) that gas would participate in res- 
onant scattering alongside stars. 

We distinguish two drivers of radial migration: when 
the angular momentum of a star is changed, whether by 
scattering at an orbital resonance or by non-resonant scat- 
tering by a molecular cloud, the star's guiding-centre radius 
changes and the star's entire orbit moves inwards or out- 
wards depending on whether angular momentum is lost or 
gained. When a scattering event increases a star's epicy- 
cle amplitude without changing its angular momentum, the 
star contributes to the density over a wider range of radii. 
In a slight modification o f the terminology introduced by 
ISellwood fc BinnevI l|2002l ). we say that changes in angular 
momentum cause "churning" while changes in epicycle am- 
plitude lead to "blurring". This paper extends models of 
Galactic chemical evolution to include the effects of churn- 
ing and blurring. 

Given the strength of the arguments that cold gas 
should participate in churning alongside stars, and that 
shocks induced by spiral structure cause gas to drift inwards, 
it is mandatory simultaneously to extend traditional chem- 
ical evol ution models t o incl ude radial fiows of gas within 
the disc. iLacev fc Fall (j 19851 ) studied chemical evolution in 
the presence of a radial inflow of gas and demonstrated that 
a radial flow enhances the metallicity gradient within the 
disc. This enhancement plays an important role in our mod- 
els, which differ from those of Lacey & Fall in that they 
include both radial gas ffows and radial migration of stars. 
Moreover we can fit our models to o bservational data th at 
is much richer than that available to iLacev fc FalU (|l985l ). 

Our models are complementary to ab-initio mod- 
els of galaxy forrnation such as those pr e sented by 
ISamland fc GerhardI (|2003l ) and iRoskar et al] l|2008bl ) in 
that they allow the solar neighbourhood to be resolved in 
greater detail, and because they are enormously less costly 
numerically, they permit parameter searches to be made that 
are not feasible with ab-initio models. 

The paper is organized as follows. Section[2]presents the 
equations upon which the models are based. These consist 
of the rules that determine the rate of infall of fresh gas, 
the rate of star formation, details of the stellar evolution 
tracks and chemical yields that we have used and descrip- 
tions of how churning and blurring are implemented. Section 



[3] describes in some detail a "standard" model of the evo- 
lution of the Galactic disc. This covers its global properties 
but focuses on what would be seen in a survey of the solar 
neighbourhood. Section|4]presents the details of the selection 
function that is required to mimic the Geneva-Copenhagen 
sample (GCS) o f solar -neigh bourhood stars p ublish ed by 
iNordstrom et al.l (|2004l 'l and iHolmberg et al.l (|2007l '). and 
explains how this sample has been used to constrain the 
model's parameters. Section [S] explains how the observable 
properties of the model depend on its parameters. Section[6] 
discusses the relation of the present models to earlier ones, 
and discusses the extent to which it is c onsistent with th e 
analysis of solar-neighbourhood data by iHavwoodl (|2008h . 
Section [7] sums up. 



2 GOVERNING EQUATIONS 

The simulation is advanced by a series of discrete timesteps 
of duration 30 Myr. 

The disc is divided into 80 annuli of width 0.25 kpc 
and central radii that range from 0.125 kpc to 19.875 kpc. 
In each annulus there is both "cold" (~ 30 K) and "warm" 
( > 10** K) gas with specified abundances {Y, Z) of helium 
and heavy elements. The "warm gas" is not available for 
star formation and should be understood to include both 
inter-cloud gas within the plane and extraplanar gas, which 
probably contains a significant fraction of the Galaxy's ISM. 
Indeed, in NGC891, a galaxy si milar to the Milky Way , of 



Indeed, m lN(_r0 89i, a galaxy si milar to tne Mnky Way , ol 
order a third of HI is extraplanar (lOosterloo et al. |2007| ) . In 
the Milky Way this gas would constitute the "intermediate- 
velocity clouds" that arc observed at high and intermediate 
Galactic latitudes (Kalbcrla & Dedcs 2008). 

Within the heavy elements we keep track of the abun- 
dances of O, C, Mg, Si, Ca and Fe. Each annulus has a 
stellar population for each elapsed timestep, and this popu- 
lation inherits the abundances Y, Z, etc., of the local cold 
gas. At each stellar mass, the stellar lifetime is determined 
by the initial abundances, and at each age we know the lu- 
minosity and colours of such of its stars that are not yet 
dead. Each stellar population is at all times associated with 
the annulus of its birth; the migration of stars is taken into 
account as described below only when returning matter to 
the ISM or constructing an observational sample of stars. 



2.1 Metallicity scale 

The whole field of chemical modelling has been thrown 
into turmoil by the discovery that three-dimensional, non- 
equilibrium models of the solar atmosphere require the 
metal ab undance of the Sun to be Zq = 0.012 — 0.014 
([Greyesse et al. 200^ rather than the traditional value ~ 
0.019. This work suggests that the entire metallicity scale 
needs to be thoroughly reviewed: if the Sun's metallicity 
has to be revised downwards, then so will the metallicities of 
most nearby stars. Crucially there is the possibility that val- 
ues for the metallicity of the ISM require revision: some val- 
ues derive from measurements of the metallicities of short- 
lived staxs_su£li_as_Bsta^ and require downward revision 
(e.g. iDafion fc Cunhall2004l ), while others are inferred from 
measurements of the strengths of interstellar emission lines. 
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and are not evidently affected by changes in stellar metallic- 
ities. If the metallicity scale of stars were lowered while that 
of the ISM remained substantially unaltered, it would be ex- 
ceedingly hard to construct a viable model of the chemical 
evolution of the solar neighbourhood. Moreover, both the 
stellar catalogue and most of the measurements of interstel- 
lar abundances with which we wish to compare our models 
are on the old metallicity scale, and unphysical anomalies 
will become rife as soon as one mixes values on the old scale 
with ones on the new. Therefore for consistency we use the 
old solar abundance Zq = 0.019 and exclude from consider- 
able metallicity values that are on the new scale. 



2.2 Star- formation law 

Stars form according to the iKennicuttI (| 19981 ) law. Specifi- 
cally, with the surface density of cold gas Eg measured in 
Mq pc~^ and t in Myr, star formation increases the stellar 
surface density at a rate 



dE. 
dt 



1.2 X 10" 



for Eg > Ecrit 
otherwise, 



(1) 



where the threshold for star formation. Edit is a parameter 
of the model and C = E"!'^® ensures that the star- formation 
rate is a continuous function of surface density. The nor- 
malisation in equation ((T} was chosen to yield the observed 
surface densities of gas and stars near the Sun. 

The stars are assumed to be distributed in initial mass 
over the range (0.1, 100) Mq according to the Salpeter func- 
tion, dA''/dM oc M~^'^^ . The luminosities, effective temper- 
atures, colours and lifetimes of these stars are taken by linear 
interpolat ion in (Y, Z) from the values given in the BASTI 
database l|Cassisi et al.ll200d) . 



2.3 Return of metals 

The nucleosynthetic yields of individual metals are in 
many cases still su bject to significant uncertainties (e.g. 
iThomas et al.lll998h : in fact models of the chemical evolu- 
tion of the solar neighbourhood have been used to constrain 
these yields (Francois ct al. 2004). 

For initial masses in the ranges 5 — 11 M© a nd 35 — 
100 M p) values of X, Y, Z,C and O were taken from iMaedeJ 
j 19921 ) using a non-linear interpolation scheme: the paper 
gives yields Flz for a low metallicity {Z = 10~*) and 
yields Yhz for a high metallicity {Z — 0.02). Guided by 
th e metallicity-dependen ce of the sizes of CO cores reported 
bv lPortinari et al.1 (|l998[ ) we take 



Y{Z) = {l-a)Yi^z + aYHz, 



(2) 



where 







for Z < 0.005 
for 0.005 < Z < 0.0075 



^ ^ 320(Z - 0.005) ... ^ ^ ..WW , . , . 

0.8 + 16{Z - 0.0075) for 0.0075 <Z < 0.02 ^ ' 
1 for Z > 0.02 

The yields of elements other than X, Y, Z,C and O from 
stars with masses in this range wer e take n from the 07?- 
FEO database of iLimongi fc Chieflll (|2008l ) with the mass 
cut set such that 0.05 M© of ^^Ni is produced; this rela- 
tively low mass cut repr oduces the Ca/F e ratio measured in 
very metal-poor stars bv lLai et all l|2008h . Stars less massive 



than 10 Mq were assumed to produce no elements heavier 
than O. For stars with mass es < 5 Mg , the yields were taken 
by linear interpolation from lMarigol (|200ll ). 

For initial stellar masses in the range 11 — 35 M© we 
used the metallicity-depe ndent yields of heavy elements from 
IChieffi fc Limon"gil (|2004l ) by linear interpolation on mass 
and metallic ity, extrapolating up to a = 1.5 or Z = 0.03 
respectively. IChieffi fc Limongil l|2004l ) used a relatively high 
mass cut, which produced 0.1 M© of ^®Ni. With our inter- 
polation the average amount of ^®Ni produced is well within 
the expected range. 

A fraction /eject of the gas ejected by dying stars leaves 
the G alaxy; we tested models with < /eject < 0.05 l|Pagell 
Il997h . Increasing /eject has the effect of reducing the fi- 
nal metallicity of the disc; in fact there is almost com- 
plete degeneracy between the values of /eject and nucle- 
osynthetic yields. In view of the evidence that star for- 
mation near the Galactic centr e drives a Galactic wind 
(|Bland-Hawthorn fc CohenI l2003l ). we set /eject = 0.15 at 
R < 3.5 kpc in models that use the accretion law ^ below. 
At all other radii we set /eject ~ 0.04. 

A fraction /direct of the ejecta goes straight to the cold 
gas reservoir of the local annulus, and the balance goes to 
the annulus's warm-gas reservoir. Setting /direct to values 
~ 0.2 has a significant impact on the number of extremely 
metal-poor stars predicted near the Sun. However, such large 
values of /direct are not well motivated physically, and in the 
models presented here /direct = 0.01 has a negligible value. 

In each timestep St a fraction St /tcooi of the "warm" 
gas (which includes extraplanar gas) transfers to the cold- 
gas reservoir from which stars form. The parameter tcooi 
is determined by the dynamics of extraplanar gas and the 
balance between radiative cooling and shock heating within 
the plane. Consequently, its value cannot be determined 
a priori from atomic physics. Increasing tcooi increases the 
mass of "warm" gas and delays the incorporation of freshly 
made metals into new stars, so the number of very metal- 
poor stars formed increases with tcooi. Although of order a 
quarter of the neutral hydrogen of NGC 891 is extraplanar 
(|Oosterloo et al.ll2007l ) , some of this gas will be formerly cold 
interstellar gas that has been shock accelerated by stellar 
ejecta. We do not model shock heating of cold gas, and re- 
plenish the warm-gas reservoir exclusively with stellar ejecta 
(from stars of every mass). Hence the warm-gas reservoir 
should be less massive than the sum of the extraplanar and 
warm in-plane bodies of gas in a galaxy like NGC 891. We 
have worked with values tcooi <; 1 Gyr that yield warm-gas 
fractions of order 10 percent. The results of the models are 
not sensitive to the value of tcooi. 

There is abundant evidence that pristine intergalactic 
gas disappeared from the intergalactic medium (IGM) long 
ago: quasar absorption-line st udies reveal an early build up 
of heavy elements in the IGM (|Pettini et al.ll2003l ). While it 
is clear that the disc formed from material that had been 
enriched by pregalactic and halo stars, it is unclear what 
abundances this material had. We take the chemical com- 
position of the pre-enriched gas to be that of the "warm" 
ISM after two timesteps, starting with 5 x 10* Mq of pris- 
tine gas. In each of the following four timesteps, a further 
1.25 X 10* M0 of gas with this metallicity is added to the 
disc. The surface density of the added gas is proportional to 
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Figure 1. The rate of mass injection by SNIa in the standard 
model (sohd black line) versus time. The broken green line gives 
the star formation rate in the solar annulus. 



of our models have sound foundations, they can usefully con- 
strain the nature of infall. 

In principle the rate and radial distribu tion of infall 
i s dete rmined by cosmology. For example, iNaab fc Ostrikeil 
(120061 ) infer it by assuming that the d isc scale length grows 
in parallel with the cosmic scale, while IColavitti et all (|2008l ) 
derive the global rate from N-body simulations. At this stage 
we feel that cosmological simulations are beset by too many 
uncertainties to deliver even a secure global infall rate, never 
mind the radial distribution of infall. In particular the extent 
of angular-momentum exchange between baryons and dark 
matter is controversial, as are the extent to which gas is 
accreted from cold infall rather than a hot corona. Moreover, 
nothing is known with any confidence about the dynamics 
of the corona. 

For want of clear inputs from cosmology we have sought 
a flexible parametrisation of infall. First we parametrise the 
global infall rate, and then the radial distribution of infall. 



(10- g(«-19'8kpc)/11.8kpc^g-fl/4kpc 

Thus the surface density is exponential with scalelength 
4kpc inside ~ Ro but tapered to zero at the outer edge 
of the grid. The existence at the outset of a warm, pre- 
enriched component of the ISM is physically well motivated 
and proves the most effective way of producing the right 
number of metal-poor stars. 

Type la supernovae are included by assuming that 7.5 
per cent of the mass in white dwarfs formed by stars of 
initial mass 3.2 to 8.5 M© ultimately explodes in type la 
supernov ae. The yields were ta ken to be those of the W70 
model in llwamoto et all (|l999l ). It is believed that the pro- 

J enitors of type la s upernovae have lifetimes of order a Gyr 
Forster et ai1l2006l ) and we have taken the mass Mwd of 
the population that survives to time t from white-dwarf for- 
mation to satisfy 

dMwD _ f for < t < 0.15 Gyr 

dt ^ I -A/wd/1.5 Gyr for T > 0.15 Gyr. ^ ' 

The rate of type la SNe is constrained by the requirements 
that (i) [0/Fe] has to fall from ~ 0.6 for the oldest stars to 
around —0.1, and (ii) [Ca/Fe] should go from about ~ 0.3 
to ^ 0. The full curve in Fig. [T] shows for the best-fitting 
model the mass-return rate as a function of time, while the 
green dashed curve shows the SFR at the solar radius. 

2.4 Inflow 

It is generally agreed that viable models of galactic chemical 
evolution require the disc to be constantly fed with gas from 
intergalactic space; inflow resolves several serious problems, 
including (i) the appearance of too many low-mctallicity 
stars near the Sun (the "G-dwarf problem", e.g. Pagel 1997), 
(ii) excessive metallicity of the current ISM, (iii) an unre- 
alistically low abun dance of deuterium in the current ISM 
(|Linskv et al.|[200^ ). Moreover, both the short timescale for 
the current ISM to be consumed by star formation and di- 
rect manifestations of infalling gas (Sancisi et al. 2008) ar- 
gue strongly for the existence of infall. Unfortunately, many 
aspects of infall are extremely uncertain. We find that the 
predictions of our models depend sensitively on how these 
uncertainties are resolved, so to the extent that other aspects 



2.4-1 Infall rate 

We have investigated two approaches to the determination 
of the infall rate. The first starts with a quantity of gas 
(8 X 10^ M0) and feeds gas into it at a rate 

M=^e-'/''i+4^e-'/''^. (6) 
Ol 02 

Here 61 ~ 0.3 Gyr is a short timescale that ensures that the 
star-formation rate peaks early on, while 62 — 14 Gyr is a 
long timescale associated with sustained star formation in 
the thin disc. We adopt Mi ~ 4.5 x 10^ M© and choose M2 
such that after 12 Gyr the second exponential has delivered 
2.6 X 10^" Mq. 

In an alternative scheme, the gas mass within the disc 
is determined a priori and infall is assumed to be available 
to maintain the gas mass at its prescribed level. We have 
investigated schemes in which the gas mass declines expo- 
nentially with time, but focused on models in which it is held 
constant at 8.4 x 10^ M©; models in which the gas mass de- 
clines exponentially produce very similar results to models 
in which the infall rate declines exponentially. 

2.4-S Distribution of infall 

We know even less about the radial distribution of the in- 
falling gas than we do about the global infall rate. In fact 
our only constraint is that the stellar disc has an approxi- 
mately exponential surface density now, and was probably 
exponential at earlier times too. Besides the star-formation 
law, the structure of the stellar disc depends on both the ra- 
dial distribution of infall and gas fiows within the disc, and a 
disc that is consistent with observations will not be formed if 
either the radial infall profile or the internal gas fiow is fixed 
without regard to the other process. Consequently, the re- 
quirement that only observationally acceptable discs be pro- 
duced requires one to develop a parametrisation that couples 
infall and flow in a possibly unphysical way. The scheme we 
have developed involves such an unphysical coupling - this 
is the price one pays for a scheme that allows one to explore 
as economically but fully as possible a range of infall profiles 
and internal flows that are consistent with the known radial 
structure of the disc. 
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We start from the assumption that the surface-density 
of gas is at all times exponential, Sg(-R) oc e"^'^^'', where 
Rd = 3.5 kpc is chosen such that with the star-formation 
law adopted above, the inner stellar disc acquires a scale 
length Rt, = Rd/ IA = 2.5 kpc sim i lar to that deter mined 
from star counts (|Robin et al.ll2003l : [juric et al.ll2008l ). Our 
value for the scale length of the gas disc is in good agree - 
ment with the value measured bv lKalberla fc Dedej (|2008h : 
3.75 kpc. Notice that we assume not only that the stellar 
disc is exponential, but that its scale length is unchanging. 
Hence we are assuming that the disc forms simultaneously 
at all radii, rather than "inside-out" . The remarkably large 
age estimated for the solar neig hbourhood (at R = 3R, ) sug- 
gests simultaneous formation (|Aumer fc Binne"vll2008l . and 
references therein) . However, our work could be readily gen- 
eralised to inside-out growth by making Rd a specified func- 
tion of time (e.g. Naab & Ostrikcr 2006), but we reserve this 
extension for a later paper. 

Our scheme for parameterising infall and flow depends 
on two parameters, /a and /b and is easiest to explain by 
considering first the limiting cases in which one parameter 
vanishes. 

Either of the algorithms of the last subsection specifies 
what the total gas mass should be at the start of a timestep: 
this is either the prescribed constant or, when equation ^ 
is used, it is the mass in the disc at the end of the previous 
timestep plus the amount that falls in during the most recent 
timestep. Hence the mass that should be in each annulus at 
the start of a timestep follows from the assumed exponential 
profile of the gas disc. Subtracting from this the mass that 
was present after the previous timestep, we calculate the 
need, i.e. the amount of gas that has to be added, of the ith 
annulus A Mi. 

We fill annuli up with gas in sequence, starting with 
the innermost ring 0. When /b = ("Scheme A") this an- 
nulus receives /aAMq from the IGM, and grabs the balance, 
(1 — /a)AMo, from annulus 1, where /a — 0.2 is a parameter 
of the model. Annulus 1 receives /aAMi from the IGM, and 
grabs the balance of its requirement, (1 — /a)(AMi -|- AMq), 
from annulus 2. The updating of every annulus proceeds sim- 
ilarly, until the last annulus is reached, which covers its en- 
tire need from the IGM. The characteristic of this "Scheme 
A" is the development of a large flux of gas through the 
outer rings - an example is given by the full red curve in 
the upper panel of Fig. [5] This flux transports inwards met- 
als synthesised in these rings and tends to deposit them at 
intermediate radii, where the inward flux is diminishing. 

When /a = ("Scheme B") annulus obtains /bAMq 
from the IGM and the rest from annulus 1. Annulus 1 now 
obtains /b[AMi + (1 - /b)AMo] from the IGM, and so on 
to the outermost ring, which is again entirely fed by the 
IGM. The short dashed blue curve in the upper panel of 
Fig. [2] shows a typical example of a mass flow through the 
disc with Scheme B. Whereas the flow generated by Scheme 
A (red curve) increases monotonically from the centre, the 
Scheme-B flow rises quickly with galactocentric distance R 
near the centre but then peaks at i? ~ 5 kpc. In the outer 
region in which the inflow is small, the metallicity forms a 
plateau. The extent of this plateau is controlled by /b: the 
larger /b, the smaller the radius at which the inflow rate 
peaks and the further in the metallicity plateau extends. 

In either of these schemes a flxed fraction of each an- 
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Figure 2. Upper panel: The rate of flow of gas over the circle of 
radius R induced by infall Scheme A with /a = 0.6 (red curve), 
infall Scheme B with /b = 0.05 (blue short dashed curve), and 
infall Scheme AB with /a = 0.35 and /b = 0.025 (green long 
dashed curve; the standard model). Lower panel: the correspond- 
ing rates of accretion from the IGM per unit area of the disc. 



nulus's need is taken from the IGM, but the definition of 
"need" is different in the two schemes: in Scheme B it in- 
cludes the gas that was taken from it by its inner neigh- 
bour, and in Scheme A it does not. In Scheme A only a 
fixed fraction of the local need is provided by the IGM, so 
the flow Fr in the disc continuously builds up through the 
disc. In Scheme B, by contrast, a part of the flow required 
in Scheme A is met by additional accretion. Consequently, 
if one wrote an equation for dJV/dr, a term — (/sAr)i<V 
would appear, where Ar = 0.25 kpc is the width of annuli, 
and this term drives exponential decay of Fr- Scheme A en- 
hances the metallicity of the middle section of the disc and 
causes the metallicity gradient to be steepest towards the 
outside of the disc. Scheme B enhances the metallicity of 
the inner disc and flattens the gradient at large radii. 

In Scheme A, if /a is set too low, the flux of gas through 
the outer annuli becomes implausibly large in relation to the 
mass of gas that is in these annuli, and radial flow veloci- 
ties ?;_B>20kms~^ are predicted. In Scheme B /b can be 
quite small because, although the flow of gas through the 
disc builds up more quickly at small radii, it peaks at a few 
kiloparsecs and then declines to small values in the outer 
disc. If either /a or /b is large, the flow through the disc be- 
comes small and the metallicity of the solar neighbourhood 
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becomes unrealistically large through the accumulation of 
metals created at the solar radius and beyond. 

Satisfactory fits to the data can be obtained only when 
both /a and /b are non-zero In this "Scheme AB" annulus 
receives a mass (/a + /b)AMo from the IGM and grabs the 
balance Men = (1 — /a — /b)AMo from annulus 1. Annulus 1 
receives a mass /aAMi +/B(AAfi-(-Moi) from the IGM and 
grabs the balance of its requirement from annulus 2, and so 
on. Notice that the radial flow profile in Scheme AB is not 
simply the sum of the corresponding profiles for Schemes 
A and B used alone. The green curve in the upper panel of 
Fig.[2]shows the radial fiow profile obtained with Scheme AB 
with the parameters of the standard model. In this model 
the radial velocity of disc gas currently rises roughly linearly 
from zero at the centre to 1.3 kms~^ at the Sun. Beyond the 
Sun a plot of radial velocity versus radius gradually steepens 
to reach 5 kms~^ at the edge of the disc. 

For each accretion scheme, the lower panel of Fig. [5] 
shows the corresponding radial distribution of accretion 
from the IGM. 



2.4.3 MetalUcity of the IGM 

We have to prescribe the metallicity and alpha-enhancement 
of gas taken from the IGM. It is far from clear how this 
should be done. 

Quasar absorption line-stud ies reveal an early build up 
of heavy elements in the IGM (jPettini et all 120031 ). More- 
over, the handful of high-velocity clouds for which metallici- 
ties have been meas ured, have heavy-element abun dances of 
order a tenth solar (Ivan Woerden fc Wakkeij[200i ). Finally, 
the metallicities of the most metal-poor thick-disc stars are 
similar to the metallicities of the most metal-rich halo stars, 
which suggests that the early disc was pre-enriched by pre- 
galactic and halo stars. We assume that throughout the sim- 
ulation accreted gas has metallicity Z = O.IZ0. 

Given that the thick disc is alpha-enhanced l|Venn et al.l 

I2OO4I '). it is clear that when disc formation starts, infalling 
gas must be alpha-enhanced. It is natural that this enhance- 
ment should decline with time as Fe from type la SNe finds 
its way into the IGM. In deed, in addition to gas that fiow s 
out in the Galactic wind (|Bland-Hawthorn &: Coh"enll2003h . 
type la SNe in dwarf spheroidal galaxies will have con- 
tributed their Fe to the local IGM, and if the Magellanic 
Stream is made of gas torn from the SMC, it will have been 
enriched with Fe from SNe in the SMC. Thus we expect the 
metallicity and alpha enhancement of the IGM to be time 
dependent and governed by the chemical-evolution histories 
of galaxies. 

These considerations suggest making the a- 
enhancement of the IGM refiect that of an outer annulus 
of the Galaxy; the chemical evolution of this ring acts 
as a proxy for the combined chemical evolution of the 
many contributors to the chemical evolution of the IGM. 
If the IGM were assumed to mirror the outermost ring, its 
Q-enhancement would remain extremely low because this 
annulus takes all its gas from the IGM, and passes what 
few heavy elements it synthesises inwards. Hence the IGM 
must mirror an outer annulus but not the outermost. In our 
models the a-enhancement of the IGM mirrors the annulus 
with radius R = 12.125 kpc. Since yields of a elements 



decline with increasing metallicity, the outer disc should be 
a-enhanced. 

2.5 Churning 

Transient spiral arms cause both stars and gas to be ex- 
changed between annuli in the vicinity of the corotation res- 
onance. Such exchanges automatically conserve both angu- 
lar momentum and mass. Since these exchanges are driven 
by spiral structure, in which hot and extraplanar gas is not 
expected to participate, churning is confined to stars and 
cold gas. We restrict exchanges to adjacent rings but allow 
two exchanges per timestep, so within a timestep second- 
nearest neighbouring rings exchange mass. 

Further studies of spiral structure in high-quality N- 
body simulations are required to determine how the proba- 
bility of a star migrating varies across the disc. In the ab- 
sence of such studies the following dimensional argument 
suggests what the answer might be. Consider the probabil- 
ity Pex that in a characteristic dynamical time (where 
K is the local epicycle frequency) a star is involved in a 
resonant exchange across corotation. It is natural that a 
process dependent on gravitational self-energy in the disc 
should scale with the square of the surface density. Toomre's 
Q = (tk/ttGE, where a is the radial velocity dispersion, is 
a dimensionless variable, so we conjecture that Pcx oc 1/Q^. 
Our grid is uniform in R whereas an exchange across coro- 
tation changes R by of order the most unstable wavelength 
Adit = gQ/ K. The number of swaps between rings required 
to wander a distance Acrit scales as A^j.it. Moreover, the num- 
ber of ring-swaps in time scales as so the ring-swap 
probability per timestep Pring should be Acrit^ times Pex- 
This argument yields Pring oc ct^/k. In realistic cases cr'^ oc E 
and K oc P, so Pring oc SP oc M, the mass of a ring. This 
argument suggests that we take the probability pij that in 
a given half-timestep a star or gas cloud in the ith annulus 
is transferred to the jth annulus to be 

^ f fcchMj /A/max for j = i ± 1 , . 

^"^^lO otherwise, ^' 

where Mj is the mass in cold gas and stars in the jth annulus 
and Mmax = maxj(A/j). This rule ensures that the mass 
transferring outwards from the ith annulus is proportional 
to AfiAfi+i, and an equal mass transfers inwards, ensuring 
that angular momentum is conserved. The constant fcch is 
the largest transition probability for any annulus in a given 
timestep. It is treated as a free parameter to be fitted to the 
data. 

The procedure for distributing the metals released by 
a population of stars born in annulus i is as follows. The 
probability that a star born in annulus i at timestep m is 
found to be in annulus j at timestep n is equal to the ijth 
element of the product matrix p^Pm+i x ■ • • pn. In practice 
we recompute p only each five timesteps and approximate 
Pm X • • • pm+4 by p^. Fig. |3] shows the extent to which the 
guiding centres of stars are changed over the lifetime of the 
Galaxy when fcch = 0.25. 

2.6 Blurring 

In addition to changing their guiding-centre radii through 
the churning process, stars oscillate around their guiding 
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Figure 3. The radial distribution of the guiding centres of 12 Gyr- 
old stars that were born at 5 (red), 7.6 (green, long dashed) and 
lOkpc (blue, short dashed) when the churning fraction k^i^ = 
0.25. 



centres with steadily increasing amplitudes. Consequently, 
stars spend time away from their guiding-centre radii. For 
simplicity, we assume in this section that the circular speed 
Vc is independent of radius and that the vertical motion can 
be ignored because it decouples from motion in the plane. 

The fraction of its time that the orbit with energy and 
angular momentum E, L spends in a radial interval {R, R + 
dR) is 



dp : 



dt 

T 



dR 



^2{E - $eff ) 



(8) 



where T = h/Q.r is the half period and ^cs{R,L) = 
+ L^/2i?^ is the effective potential. We need to av- 
erage this over all stars with given L. These stars have some 
distribution over the energy E = + $cff • It is expedient 
to decompose E into the energy ^c{L) = $(-Rc) + Li^ /2Rl 
of the circular orbit (with radius Rc) of angular momentum 
L and the random energy £ = E — ^c- Following IShul (Il969l ') 
we take the distribution function (DF) to be 



e/a 



(9) 



where F{L) is a function to be determined. The DF ([9]) 
ensures that the radial velocity dispersion is approximately 
(but not exactly) a. Normalizing / such that J dLdJu f = 
1, where Jr{£) is the radial action, the probability that a 
randomly chosen star lies in {R, R + dR) is J dLdJn fdp. 
Recalling that dLdJn — dL d£ /Qr and substituting for / 
and dp, we find that the number of stars in the annulus is 



dn{R) = dL djR (/dp) 



dL^ 



NdR 



NdR 



-e/a 



"I>c - "I>c 



=(10) 



da; ■ 



where N is the total number of stars in the system. The 
integral over x is simply a J dte~* /y/t = y^a. Thus we can 
conclude that the probability per unit area associated with 
a star of given L is 



PiR) 



dn 



K 



N2TvRdR aR 



exp 



ML)''i>cttiR,L) 



(11) 



where K is chosen such that 1 = 2-k J dR RP{R). 

The parameter a used in these formulae is actually 
smaller than the rms radial velocity dispersion, which is 
given by 



(4) 



7?E 



dLFaexp[($c - $off)/cr^ 



where the stellar surface density is 



R 



E 

dL — exp[($c 
a 



$cff)/a' 



(12) 



(13) 



For specified radial dependencies of E and {v%), equations 
(|12p and (I13II can be use d to determine the functions F{L) 
and (t(L) lDehnenlll999l ). However, in the present applica- 
tion it is not E(7?) that we wish to specify, but the number 
of stars with guiding centres in each ring: 



diV „ dN 

VcENtot f d£_ 

^R 



djR f{L, Jr) 



(14) 



where A^'tot is the total number of st ars in the disc. We adapt 
the technique described by Dehnen (1999) for determining 
F{L) and a{L) from equations p2|) and (|13|) to the deter- 
mination of these quantities from equations (|12p and (|14|) . 
Specifically, we start from the values of F{L) and cr{L) that 
would hold in the epicycle approximation, when i}R = k 
independent of £ and 



F{L) 



dN 



iiciVtot diic 



(15) 



Then at each L we evaluate (ii^) from (|12p and multiply 
a by the ratio of the desired value to the value just cal- 
culated. Then we re-evaluate F from (|14p and repeat until 
convergence is obtained. 

We now address the question of how {vr) should depend 
on radius. The scale heights h of galactic discs are found 
to be largely independent of radius I van der Kruit fc Searld 
Il982l ). and for h R (when the vertical dynamics can 
be considered one-dimensional) this finding implies that 
the vertical velocity dispersion scales with the surface den- 
sity as Y}^^ . If the ratio of the vertical and radial ve- 
locity dispersions (Jz/i yR)^^^ is independent of rad ius, as 
is often assumed (e.g. iKregel fc van der KruitI |2005^ . then 
(i;^) oc S oc e~^^^* , In the solar neighbourhood at i? ~ 3R* 
the oldest stars have cr_R > 40 km s~^, so this line of reason- 
ing predicts that {vr)^^^ > 180kms~^ in the central regions, 
which is implausibly large. 

Evidently these naive arguments based on complete 
decoupling of planar and vertical motions are inadequate 
for the old disc; we need a distribution function that 
treats the third integral properly. Pending the availabil- 
ity of such a DF we have adopted the assumption that 
(i;^) oc e~^/i'^^* , which implies that at Rt the old disc has 
(^fl)^''^ — 85kms~^, which is only sligh tly lower than th e 
velocity d ispersion in the Gala ctic bulge (iRich et al.ll2007l ). 

From lBinnev et al. I (|2000| ) we take the time dependence 
of {v^) (Ro) at solar galactocentric distance Ro to be 
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Figure 4. The radial distributions of stars with guiding centres at 
5 (red), 7.6 (green) and lOkpc (blue) when at the Sun = 
25kms~^ (full curves) and 40kms~^ (dashed curves). 



^(v%){Ro,t) = max \ 10, 38 



t + 0.038 Gyr 
10.038 Gyr 



kms 



(Tz(r) = max < 4, 25 



10 Gyr 



kms 



Further assuming that stars of a given age form an isother- 
mal population, their vertical density profile is 



n{z) 



-i,{z)/o 




which is consistent with the data of lHolmberg et ahl l|2007h . 

Fig.|4]shows blurring distributions P{R) from equation 
(jlip for three radii (5, 7.6 and lOkpc) and two values of 
(v'ji)^^^ at the Sun, namely 25kms~^ and 40kms~^. 

Note that scatterings by spiral arms and molecular 
clouds that heat the disc, also change the angular momen- 
tum of each star and therefore its guiding centre. Hence such 
scatterings contribute to both churning and blurring. 

Since churning moves the guiding centres of the stars 
themselves we first apply the churning matrix and apply 
the blurring matrix afterwards. 



2.7 Vertical structure 

For comparison with observations of the solar neighbour- 
hood we need to know the vertical distribution of stars near 
the Sun. We determine this by adopting a r elationship be- 
tween time and vertical velocity dispersion (|Binnev et al. I 
I2OOOI ) 



Figure 5. The metallicities of the current ISM in the stan- 
dard model, [Z /H] (red dashed curve) and [O/H] (red full 
curve), as functions of Galactocentric radius. Also measure- 
ments of the metallicities of HII regions by Shaver et al. (1983) 
(dark blue), Vilchez & Esteban (1996) (black crosses) and RoUe- 
ston et al. (2000) (light blue crosses). The green line shows 
(16) the linear least-squares fit to the measurements: it has a slope 
of —0.082 dex kpc~^ . The data points have been updated and 
rescaled to Rq = 7.5 kpc as described in the text. Where neces- 
sary points have been shifted vertically by —8.93 to put them on 
the solar scale. 



(17) 



(18) 



where $(z) is the difference in the gravitational potential 
between height z and the plane. This pot ential is calculated 
for a m odel similar to those presented bv lDehnen fc BinnevI 
(|l998h but with the thin and thick disc scaleheights taken 
to be 0.3 and 0.9 kpc, the total stellar surface density set to 
35.5 M© pc~^ with 3/4 of the stellar mass in the thin disc, 
and the g as surface de n sity s et to 13.2 Mp) pc~ ' ^ in c onfor- 
mity with iFlvnn et all l|2006l ) and ljuric et al.1 (I2OO8I). The 
disc scalelength is taken to be Rd = 2.5 kpc (|Robin et al.l 
|2003| ) and the dark halo density is set such that Vc{Ro) = 
220 kms"\ 



3 THE STANDARD MODEL 

We now describe the properties of our standard model as a 
preliminary to explaining how these properties depend on 
the input assumptions and the values of the various param- 
eters. In the standard model the accretion rate is given by 
equation (|6]); the values of the parameters for this model are 
given in Table [T] 

The red dashed curve in Fig. [S] shows the current 
metallicity Z of the ISM as a function of radius. There is 
quite a steep outward decline in metallicity, the gradient 
in the vicinity of the Sun being of order — 0.11 dexkpc"^. 
The solid red curve shows that [O/H] falls less steeply 
with R than does [Z/H], having a gradient near the Sun 
~ —0.083 dexkpc"^. The shallower gradient in oxygen re- 
flects our use of metallicity-dependent y ields. Although 
shallower gradients are generally cited (e.g iRoUeston et al.l 
|200Q) the data points in the figure are con s istent with the 
model. The data derive from IShaver et al.1 ll 19831) who as- 
sum ed Rp = 10 kpc and from IVilchez fc EstebanI (|l996l ) 
and iRoUeston et al. 1 l|2000l ). who assumed Ro = 8.5 kpc. To 
plot these data on a consistent scale with _Ro ~ 7.5 kpc 
we have when possible recalculated the Galactocentric dis- 
tances from the he li ocentr ic dis tances, taking the latt e r from 
iKharchenko et al.l (|2005l ) or iLoktin fc Beshenov ' (200± 
when possib le. For some of t h e poin ts in Vilchez fc EstebanI 
l| 19961 ) and iRolleston etlll (|2000l ) heliocentric distances 
were not available, so we simply reduced the cited Galacto- 
centric distance by 1 kpc. The green line in Fig.[S]is the linear 
least-squares fit to the data; its slope is — 0.082 dexkpc"^. 
Our gradient in [O/H] lies within the range of frequently oc - 
currmg values in Table 4 of lVila-Costas fc Edmund^ l|l992l ). 
who assembled data for 30 disc galaxies. 

The upper panel in Fig. [6] shows the evolution of Z 
for the cold ISM in a number of annuli - the solar annu- 
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± CXI CilllU, Lt/J. 


A/T P n Ti 1 VI o" 






^crit 


Kennicutt's threshold surface density 


limited impact 





Mo 


initial gas mass 


affects only N{Z) at [Z < —0.7; from Hess diagram 


3.0 X 10^ Mq 


Ml 


early infall mass 


affects only N{Z) at [Z < —0.7; from Hess diagram 


4.5 X lO'^ Mq 


Ma 


long timescale infall mass to 12 Gyr 


fixed by present mass 


2.9 X 10^0 Mq 


bi 


early infall timescale 


affects only N{Z) at [2/H] < —0.7; from Hess diagram 


0.3 Gyr 


62 


long infall timescale 


limited impact; estimated from Hess diagram and other work 


14 Gyr 


/a 


Scheme A fraction of gas from IGM 


free parameter 


0.36 


/b 


Scheme B fraction of gas from IGM 


free parameter fixed by local gradient 


0.025 


^ch 


churning amplitude 


free parameter 


0.35 


to 


delay before first type la SNe 


taken from literature 


0.15 Gyr 




timescale for decay of type la SNe 


taken from literature 


1.5 Gyr 


/cj cct 


fraction of ejecta lost to Galaxy 


small impact; mainly affects metallicity scale 


0.15 - 0.04 


/direct 


fraction of ejecta to cold ISM 


small impact limited to [Z /H] < —0.7 


0.01 


^cool 


cooling time of warm gas 


fixed by present mass of warm gas 


1.2 Gyr 




initial warm gas mass 


impact limited to N{Z) at [Z /R] < -0.7 


5 X 10* Mq 


ZlGM 


mctallicity of the IGM 


limited to H > 12kpc; taken from literature 


O.IZq 



Table 1. Parameters of the standard model. The infall rate is given by equation The larger value of /eject applies at i? < 3.5 kpc. 
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Figure 6. Upper panel: the metallicity of the cold ISM in each 
annulus as a function of lookback time showing each fifth ring. 
The curve for the solar annulus is red. Lower panel: the present 
density of solar-neighbourhood stars in the age-metallicity dia- 
gram. The colours encode the logarithm of the density of stars. 



lus is coloured red. The smaller the radius of an annulus, 
the higher its curve lies in this plot because chemical evo- 
lution proceeds fastest and furthest at small radii. At small 
radii the metallicity of the cold ISM continues to increase 
throughout the life of the Galaxy, whereas at _R > Rq kpc. 



Z peaks at a time that moves earlier and earlier as one 
moves out, and declines briefly before flattening out. This 
phenomenon reflects a combination of dilution by infalling 
metal-poor gas and the inward advection of metals by the 
flow through the disc. 

The lower panel in Fig. [6] shows the corresponding 
present-day metallicity distribution of solar-neighbourhood 
stars. Although this is the distribution of stars currently 
in the solar annulus, it is clearly made up of a series of 
curves, one for each annulus in the model. The curves for 
interior annuli go from green to yellow as one goes forward 
in time, reflecting the fact that relatively recently formed 
stars are much less likely to have moved a large radial dis- 
tance than older stars. Similarly, in the bottom part of the 
figure the colours go from blue to green to yellow as one 
moves towards the time axis, because then one is moving 
over curves for larger and larger radii, where both the star- 
formation rate and the probability of scattering in to the 
solar radius are low. Hence, regardless of stellar age, most 
solar-neighbourhood stars have Z in a comparatively narrow 
range centred on [Z/H] ~ —0.1. 

The full green curve in Fig. [7] shows the surface density 
of the stellar disc at 11.7 Gyr, which is roughly exponential. 
The red points show what the surface density would be if 
stars remained at their radii of birth. By construction this 
forms an exponential disc with a scalelength of 2.5 kpc. The 
broken green line shows that at R < Ro the disc approxi- 
mates an exponential with a larger scale length ~ 2.8 kpc. 
The blue curve shows the surface density contributed by 
stars formed in the first 0.8 Gyr, which will be a-enhanced. 
This distribution deviates more strongly from an exponen- 
tial because radial migration is most important for old stars. 
Fitting an exponential to this curve at R < Ro yields a scale- 
length 3.1 kpc. 

Fig. |8] reveals that thin and thick disc components can 
be identified within this overall envelope: the upper panel 
shows that the vertical stellar density profile at the Sun is 
not exponential but can be fitted by a sum of two expo- 
nentials. There is significant latitude in these fits and the 
fraction of stars that is assigned to each component varies 
wi th their sca l elengt hs. For comparison with recent results 
of lJuric et ai] (|2008h we present a fit with their value for 
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Figure 7. Full green curve: the surface density of the stellar disc 
at 11.7 Gyr. Broken green line: exponential fit to the inner part of 
this curve. Red curve: the surface density if stars remained where 
they were born. Blue curve: surface density contributed by stars 
born in the first 0.8 Gyr. Broken blue line: linear fit to this curve. 




10 12 
R/kpc 

Figure 8. Upper panel: the volume density of stars at R= 7.6 kpc 
as a function of height (red), a fit (black dashed) and its decom- 
position into thin and thick components. Lower panel: the volume 
density of stars at the current epoch at z = (red), z = 0.75 kpc 
(green dashed) and z = 1.5 kpc (blue short dashed). 
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Figure 9. The predicted distribution of solar-neighbourhood 
stars in the ([Fe/H], [O/Fe]) plane. The sample is obtained by 
using the selection function of the GCS survey as described in 
Section |4] below. The colours of points depend on the star's az- 
imuthal velocity: < 179kms~l blue; 179 < v^/kms~^ < 244 
red; > 244kms~^ green. The black curve shows the trajectory 
of the solar annulus. 



the local thick disc fraction of 13 per cent. This yielded 
scaleheights of hi — 335 pc and /12 = 853 pc, very well in 
the range of their results. Note that the double-exponential 
density structure is not caused by any pecularity in star for- 
mation history, like a peak in early star formation, but is a 
consequence radial mixing combined with the given vertical 
force field. However, precise characterisation of the vertical 
structure must await dynamical models that employ a more 
accurate form of the third integral of galaxy dynamics. 

The lower panel of Fig. |8] shows that ai z — 1.5 kpc 
(where the thick disc is dominant) the stellar distribution 
is less centrally concentrated than it is in the plane; if one 
were to fit an exponential profile to the stellar density at 
z — 1.5 kpc for R < 10 kpc, the scalelength fitted would be 
larger than that app ropriate in t he plane. Just this effect is 
evident in Fig. 16 of ljuric et al] (fioOS ). 

Fig. |9] shows the predicted distribution of solar- 
neighbourhood stars in the ([O/Fe], [Fe/H]) plane when a 
sample is assembled using the GCS selection function de- 
scribed in Section |4] below. Two ridge-lines are evident: 
at top left of the figure a population starts that stays at 
[O/Fe] ~ 0.6 until [Fe/H] ~ -0.75 and then turns down 
towards (0,0), while a second larger population starts at 
about (-0.75,0.25) and falls towards (0.2,-0.05). This ar- 
rangement of_£oints is very similar to that seen in Fig. 2 of 
IVenn et al.l (|2004l ). The upper ridge-line is associated with 
the thick disc, and the lower ridge-line with the thin disc. 
In the appendix we show that such bimodal distributions in 
[O/Fe] are a natural consequence of the standard assump- 
tions about star-formation rates and metal enrichment that 
we have made. The structure is not a product of the double- 
exponential nature of the standard model's infall law; the 
model with a constant gas mass displays exactly the same 
structure. Breaks in the Galaxy's star-form ation history 
llChiappini et al1ll997l ) and accretion events l|Bensbv et al.l 
I2OO5I I have been hypothesised to account for the dichotomy 
between the thin and thick discs. Our models reproduce the 
dichotomy without a break or other catastrophic event in 
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Figure 10. Upper panel: the distribution of GCS stars in the 
([Z/H],!)^) plane. Lower panel: the prediction of the standard 
model. Colours and contours reflect the density on a logarithmic 
scale with a 0.2 dex spacing for contours. 



our model's star-formation history. When comparing Fig. [9] 
with similar plots for observational samples, it is important 
to bear in mind differences in selection functions: Fig.[2]is for 
a kinematically unbiased sample, while most similar obser- 
vational plots are for samples that are kinematically biased 
in favour of "thick-disc" stars. 

The full curve in Fig.|5]shows the trajectory of the solar- 
neighbourhood ISM. At low [Fe/H] this runs along the ridge 
line of the thick disc, and it finishes on the ridge line of the 
thin disc, but it is distinct from both ridge lines. The sharp 
distinction between this curve and the ridge line of the thin 
disc make it very clear that the latter is formed through the 
migration of stars into the solar neighbourhood, not through 
the chemical evolution of the solar neighbourhood itself. In 
many previous studies it has been assumed that the ridge 
line of the thin disc traces the historical evolution of the local 
ISM. Fig. [9] shows that this assumption could be wrong and 
that inferences regarding the past infall and star-formation 
rates that are based on this assumption are not to be trusted. 

In Fig. [5] the points are colour coded by their an- 
gular momenta/guiding centres: blue points are for < 
179 kms^^ (i?g < 0.81i?o), red points are for 179 kms"^ < 



< 244kms~^ and green points are for > 244kms~^ 
(Rg > 1.17?o)- At the low-metallicity end of the thin-disc 
ridge line many points are green and few blue, while at the 
high-metallicity end the reverse is true. Thus low-metallicity 
thin-disc stars tend to have guiding c entres Rg > Rp , w hile 
high metallicity stars have Rg < Ro- iHavwoodl l|2008l ) has 
noted the same metallicity-velocity correlations in samples 
of nearby stars. The thick disc contains stars from all three 
radial ranges, but stars with small Rg (blue) are most promi- 
nent at higher [Fe/H]. 

Fig. llOl shows the distribution of stars in the ([Z/H], ^0) 
plane: the upper panel is for the GCS stars and the lower 
panel is for the standard model. In both panels the highest 
density of stars lies near (0, 220 kms~^) and the upper edge 
of the distribution rises as one moves to lower metallicities. 
The metallicity gradient in the disc leads to the main clus- 
ter of stars sloping downwards to the right. A significant 
difference between the two panels is that in the upper panel 
there are more stars in the lower left regio n. This populatio n 
is very much more prominent in Fig. 5 of IHavwoodl l|2008h . 
where a band of points runs from small and [Z/H] up 
towards the main cluster. This band is made up of halo and 
thick-disc stars that are selected for in the samples from 
which Haywood drew data. The other difference between 
Haywood's Fig. 5 and the lower panel of Fig. [10] is that 
Haywood's main clump has a slightly less pronounced slope 
down to the right. It is likely that errors in the measurements 
of [Fe/H] have moderated this slope. The GCS distribution 
shown in the upper panel of Fig.[TO]is (especially on the high 
metallicity side) dominated by overdensities around the ro- 
tational velocitie s of well-know n stellar streams (e.g., the 
Hercules stream lDehnen|[l998l ). This pattern overlays the 
general downwards slope. The model accounts well for the 
steeper edge of the density distribution at high rotational 
velocities, which is the combined effect of lower inwards blur- 
ring and lower stellar densities from outer rings. 

The top panel of Fig. [TT] shows that the stellar metal- 
licity distribution is less centrally concentrated than that of 
the cold ISM from which stars form. Three factors are re- 
sponsible for this result. First, the mean metallicity of stars 
reflects the metallicity of the gas at earlier times, which 
was lower. This effect is most pronounced at the centre, 
where the metallicity of the ISM saturates later than fur- 
ther out. Second, radial mixing, which flattens abundance 
gradients, has a bigger impact on stars than gas because 
stars experience both churning and blurring. Third, the net 
inflow of gas st e epens the abundance gradient in the gas. 
iHolmberg et al.l l|2007l ) have estimated the stellar metallic- 
ity gradient from the GCS stars. When they select thin- 
disc stars they find —0.09 dexkpc"^, but when one excludes 
stars with [Z/H] < —0.7 (which ensures halo objects are 
removed), one obtains —0.11 dexkpc"^. The gradient of the 
dashed red line in Fig. [11] at 7.6 kpc is O.lOdexkpc"^ in 
excellent agreement with the GCS data. 

The lower panel of Fig. [TT] shows the breadth of the 
metallicity distribution at three radii. These distributions 
have full-width at half maximum around 0.35 dex and are 
significantly offset to each other by 0.25 dex. 

Fig. [12] shows how a-enhancement varies in time and 
space, in stars and gas. Naturally, [O /Fe] declines with time 
in both the ISM and in the stellar population, and at a given 
time is higher in the stars than the gas. [0/Fe] generally 
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Figure 11. Upper panel: mean metallicities of stars (dashed) and 
cold ISM (full) as functions of R at the present time (red) and 
at 1.5 Gyrs (blue). Lower panel: the distributions over metallicity 
of stars currently at ij = 5kpc (blue), 7.6 kpc (red), and lOkpc 
(green). 
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Figure 12. Full lines: [O/Fe] in the cold ISM after 1.5 Gyr (blue) 
and 12 Gyr (red). Dashed lines: Mean [O/Fe] of stars after 1.5 Gyr 
and 12 Gyr. 
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Figure 13. Absolute magnitudes in the B (blue short dashed), 
R (red) and / (black long dashed) bands as functions of time. No 
allowance has been made for obscuration. 

increases outwards but at 12 Gyr in both stars and gas it 
attains a plateau at J?> 10 kpc, with [a/H] ~ 0.2 in the gas. 
The existence of the plateau is a consequence of the rule 
that in the IGM [0/H] is the current value in the disc at 
i? ~ 12 kpc; gas with the given a-enhancement rains on the 
disc at i?<20kpc, is enriched by supernovae of both types 
and a few gigayears later arrives at iZ = 12 kpc with its 
original a-enhancement. This level is set by the metallicity- 
dependent yields we have employed. 

Fig. [13] shows the B, R and /-band absolute magni- 
tudes of the standard model as functions of time. The B- 
band luminosity rises quickly to a shallow peak around 5 Gyr 
and then commences a very slow decline. Emissions in the 
i?-Band are almost constant at the present time, while I- 
band luminosities continue to rise throughout the Galaxy's 
life because additions to the stock of long-lived stars out- 
weigh deaths of relatively short-lived and predominantly 
blue stars. In our model the Galaxy reaches an I-Band 
magnitude of around —22.7 which is exactly the result one 
would expect fo r a disc galaxy with a rotation velocity of 
220 km s"^ (e.g. IPizaeno et ahboosh 

Fig. ll4l shows the U — B and B — I colours of the disc at 
t = 1.5, 4.5, 8.4 and 12 Gyr. As expected, the disc reddens 
at a declining rate throughout its life. There is at all times 
a significant colour gradient between R — 8 and 16 kpc that 
make the disc's edge about 0.2 mag bluer in B — J than its 
centre. 

Fig. [15] shows that radial migration causes the depen- 
dence of velocity dispersion on time for stars that are cur- 
rently in the solar neighbourhood to differ materially from 
the acceleration law (|17p that determines the time depen- 
dence of az for stars that are born at given radius. The out- 
ward migration of stars brings to the solar neighbourhood 
stars that carry with them the large velocity dispersions 
characteristic of their places of birth. The impact that these 
migrants have on for stars of a given age increases with 
age, so at high ages increases faster than equation (|17p 
predicts. Least-squares fits of oc t'^ to the red and green 
curves in Fig. 1151 vield l3 — 0.53 and 0.44, respectively. Em- 
pirically the photometrically complete portion of the Hip- 
parcos catalogue shows that the best power-law fits to the 
rate of increase of az yield /3 ~ 0.45 (|just fc Jahreissll2007l : 
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Figure 14. U — B (upper panel) and B — I (lower panel) as 
functions of radius at t = 1.5, 4.5, 8.4 and 12Gyr. No allowance 
has been made for obscuration. 



lAumer fc Bin"^l2008l ) . From a theoretical standpoint, this 
result has hitherto been puzzling because the largest expo- 
nent that ca n be obtained from the dynamics of star scat- 
tering is 1/3 jBinnev fc Lacevlll988l ). Such studies treat the 
acceleration as a local process. Our result suggests that the 
conflict between theory and observation is attributable to 
violation of this assumption. 



4 FITTING THE MODEL TO THE SOLAR 
NEIGHBOURHOOD 

A major constraint on the models is provided by comparing 
the model's predictions with samples of stars observed near 
the Sun. To make these comparisons, we have to reproduce 
the selection functions of such samples, which proves a non- 
trivial job. 

The GCS is an important sample, an d for each model 
we ca lculate the likelihood of this sample. iNordstrom et all 
l|2004 ) obtained Stromgren photometry and radial velocities 
for a magnitude-limited sample of 16 682 F and G dwarfs, 
nearly all of which have good Hipparcos parallaxes. From the 
photometry they estimated metallicities and ages. There has 
been some debate about the calibration of the metallicities 
and ages ijHavwoodl [200^ : iHolmberg et al.l 
I2OO81) . Recently the re-calibrated data from 
1120071) became available and it is to these data that 
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Figure 15. Velocity dispersion of solar-neighbourhood stars as a 
function of age. Red curve: cr^ for all stars in the solar annulus. 
Green (long dashed) curve: for stars within 100 pc of the Sun. 
Blue (short dashed) curve: cr^ for stars born in the solar annulus. 



have compared our models. We compare their metallicities 
([Me/H]) to our [Z /H] as it is not entirely obvious to what 
extent alpha enrichment enters into their measurements. 
Since assigning ages to individual stars is very difficult, we 
have concentrated on matching the distribution of stars in 
the {MvjTeB, Z) space from which ages are derived. 

For each metallicity we construct a volume-limited stel- 
lar number density of stars in the (MvjT^a) plane by con- 
sidering each annulus j, and calculating the fraction of each 
population in this annulus that will be in the solar neigh- 
bourhood. For given absolute magnitude, the probability 
that a star will enter the sample is 



drr / d Qn{z), 



(19) 



where the space density of stars n{z) is assumed to be plane 
parallel and given by equation ((18} . 

For the GCS selection funct ion we use the ap proximate 
b — y colour rules from ,Nordstrom et all (|2004| ) - a more 
sophisticated selection function could in principle be con- 
structed, but it is not possible from the published data. At 
each colour, the appropriate selection function is charac- 
terised by two apparent magnitudes vi and 112 listed in Table 
[5] <!> declines linearly from unity at magnitudes brighter than 
vi to zero fainter than V2 . (j) vanishes for b — y bluer than 
0.21. For b — y redder than 0.38, (j) is reduced by a factor 0.6. 

The selection of INordstrom et al.l (|2004l ) is designed to 
exclude red giant branch (RGB) stars from the data set. 
However, some of these stars are still in the sample, while the 
sample is biased against stars just below the giant branch. 
We take out of consideration the RGB itself and we down- 
scale the theoretically expected population density near the 
starting point of the red giant branch by a factor of 4 to 
reconcile it with the data. Since the number of RGB-stars 
in the GCS is not large anyway, the loss of information is 
small. We also removed from the dataset three objects that 
are far too faint to be attributed to the main sequence. The 
theoretical distributions are convolved with a Gaussian of 
dispersion 0.1 dex in [Z /H] to allow for measurement errors. 

Fig. [16] compares predicted (lower panel) and observed 
(upper panel) Hess diagrams for the GCS stars. As discussed 
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Table 2. Magnitudes defining the Geneva— Copenhagen selection 
function 

b-y 0.21-0.25 < 0.344 < 0.38 < 0.42 > 0.42 
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Figure 16. Comparison of the observed (upper) and predicted 
(lower) distributions of GCS stars in the {T^ft, My) plane. Con- 
tours have an equal spacing of 10 counts per bin, starting with 3. 



by iHolmberg et Zl JioOTl ). the ridge-line of the main se- 
quence in the GCS data is significantly displaced from that 
predicted by isochrones. We have eliminated the effects of 
this offset on Fig.[T5]in the simplest possible way, namely by 
decreasing all model values of log(rcff ) by the value, 0.015, 
that yields the closest agreement between the theoretical 
and observation main sequences. After this correction has 
been made, the agreement between the theoretical and ob- 
servational Hess diagrams shown in Fig. [16] is convincing 
though not perfect. The original conception had been to 
determine the model's parameters by maximising the likeli- 
hood of the GCS stars in the model density in (My , Tcff , Z) 
space, but confidence in this plan was undermined by (i) the 
need for an arbitrary alignment of measured and theoretical 
values of Tcn, and (ii) the extent to which the likelihood of 




-0.5 
[Z/H], [Me/H](3cs 

Figure 17. The metallicity distribution of GCS stars (green 
points) and the corresponding prediction of the standard model 
(full red curve). The broken blue curve shows the model that dif- 
fers from the standard model only in the elimination of churning 
and radial gas fiows. The broken pink curve shows the model with 
neither churning nor blurring. The lower panel shows the same 
data but with a logarithmic vertical scale to reveal structure in 
the wings of the distribution. 



the data depends on the uncertain GCS selection function. 
Notwithstanding these reservations, we are encouraged that 
the standard model maximises the likelihood of the data 
at an age, ~ 11 Gyr, that agree s with other estimates of the 
age of the solar neighbourhood l|Aumer fc Binnev„2008 . and 
references therein). 

The full red curve in Fig.[l7]shows the metallicity distri- 
bution of stars predicted by the standard model, and green 
points show the GCS data. The agreement is excellent. At 
the lowest metallicities theory predicts slightly too few stars, 
but the uncertainties in both the theory and the data are 
large in this limit and the theory does not include halo stars, 
so it should under-predict the data at [Z/Yi] ^ — 1. 



5 GENERAL TRENDS 

We now discuss aspects of how the observable properties of 
a model depend on its parameters, and what is required to 
achieve fits of the quality seen in Figs [S] and 1171 

Table [T] lists the model's 16 parameters; the third 
column explains whether the parameter is fitted to the 
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GCS data or taken from the literature, and indicates the 
sensitivity of models to that parameter. Five parameters 
(Mo, Ml, &i, /direct, Mwarm) significantly affect only the dis- 
tribution of stars at [Z /B] < —0.7. Six other parameters 
ooi, Zigm) are fixed by observed properties of 
the Galaxy other than the solar-neighbourhood stellar dis- 
tribution. The remaining five parameters are the critical sur- 
face density for star formation Ecrit , the long infall timescale 
62, the accretion parameters /a and /b and the churning 
strength kch- We shall see that /b is effectively set by the 
metallicity gradient in the gas, that 62 is effective deter- 
mined by the local Hess diagram, and that the value of Ecrit 
is unimportant providing it is small (we have included this 
parameter only for consistency with earlier work; the mod- 
els do not want it). Consequently the fit of the model to the 
data shown by the red curve and the green points in [17] is 
obtained by adjusting just /a and kch- 

The number of stars more metal poor than [^/H] ~ 
— 1 depends sensitively on the thermal structure of the 
early ISM. Most previous studies (exceptio ns include 
I Thomas et al. I Il998l : ISamland fc GerhardI l20oi l have used 
only one phase of the ISM. Introducing the warm component 
of the ISM delays the transfer of metals to the star-forming 
cold ISM by ~ 1 Gyr, thus increasing the number of ex- 
tremely metal-poor G dwarfs. Our first models initially had 
no warm gas, with the result that at early times the mass 
of warm gas was proportional to time and the metallicity of 
the cold gas rose quadratically with time. These models had 
an over- abundance of very metal-poor stars. These experi- 
ments led to the conclusion that pregalactic and halo stars 
endowed the disc with warm, metal-rich gas at the outset. 
Even at late times, the existence of the warm ISM delays 
the introduction of freshly-made metals into stars, and thus 
in concert with the gas fiow through the disc steepens the 
metallicity gradient in the stellar disc; eliminating the warm 
component raises the metallicity of the solar neighbourhood 
and beyond by ~ 0.1 dex. 

The metallicity of infalling gas only affects the struc- 
ture of the disc at R>12 kpc. Lowering Zigm steepens the 
metallicity gradient at i? > _Ro- 

It is instructive to consider the case in which blurring 
is included but churning is turned off by setting kch = and 
radial flows are eliminated by ensuring that /a + /b = 1 
so every ring's need is fully supplied from the IGM. The 
broken blue curve in Fig. [T7] shows the present-day metallic- 
ity distribution that this model predicts for the GCS. The 
peak of the distribution is much narrower than in the stan- 
dard model, and there is a striking deficiency of metal-rich 
stars. The broken pink curve shows the effect of also turning 
off blurring: the deficiency of metal-rich stars becomes even 
more striking but there is negligible change on the metal- 
poor side of the peak. 

Reducing the current SFR by making the infall rate a 
more rapidly declining function of time shifts the peak of the 
distribution to higher metallicities by reducing the relative 
strength of recent inflow and thus the supply of fresh, metal- 
poor gas. The use of an IMF that is steeper in the low-mass 
region, as has been suggested by some studies, reduces the 
mass of metals that is locked up in low-mass stars and in- 
creases the metal production by each generation. The metal- 
licity in such a regime is accordingly higher, but the shape of 
the distribution does not change. Since changes in the IMF 



at the high-mass end, which is equally uncertain, can pro- 
duce compensating variations in yields, loss rates, etc., we 
stayed with the traditional approach via the Salpeter IMF. 

In all interesting models the metallicity of the local ISM 
saturates early on. The saturation level depends on the pat- 
tern of gas flow through the disc, and on the current SFR 
relative to the mean rate in the past: the faster the decline 
in the SFR, the higher the current metallicity. Naturally the 
stars of the solar neighbourhood are on the average younger 
in models with a constant gas mass than in models in which 
the infall rate is declining according to equation (|6]). This 
relative youth is reflected in the structure of the local Hess 
diagram. We reject the model with constant gas mass be- 
cause it assigns a signiflcantly smaller likelihood to the Hess 
diagram of the GCS stars than does a model based on equa- 
tion ((6|. 

5.1 Selecting the standard model 

All our models have quite strong metallicity gradients in 
both stars and gas (Figs [5] and [11}. Since the metallicity of 
the central gas is enhanced by radial gas fiow, and models 
with large /b have larger central flows than models with 
large /a and vice versa at large radii (Fig. [2)1 , enhancing /b 
steepens the metallicity gradient at small R and diminishes 
it at large R. 

Eliminating churning and radial gas flows (by setting 
kch = 0, /a + ,/b ~ 1) dramatically reduces the metallic- 
ity gradient within both the stellar and gas discs: at the 
present epoch the gradient in the gas near the Sun falls from 
— 0.11 dexkpc"^ to —0.01 dex kpc~^. Increasing the churn- 
ing amplitude kch both increases the width of the peak in 
the predicted solar-neighbourhood metallicity distribution 
fFigs llll and ll7p and reduces the gradient in the mean metal- 
licity of stars at R<Ro. 

Fig. llSl shows how the likelihood of the GCS metallicity 
distribution plotted in Fig. [T7] varies with /a, /b and kch- 
In the upper panel favoured models (with large symbols) lie 
along a line that slopes down and to the right. Along this 
line a decrease by 0.01 in /b is compensated by an increase 
in /a by ~ 0.08. 

As one moves down the upper panel of Fig. 1181 the 
steepness of the metallicity gradient near the Sun increases, 
and the models with /b — 0.015 have local gradients steeper 
than — 0.12 dexkpc^^, which may confiict with the data. 
Models higher up the panel have smaller local metallicity 
gradients and require larger values of fcch to bring a sufficient 
variety of stars to the solar neighbourhood. Models to the 
right of the panel have smaller inward flows of gas, leading 
to local metallicities that rise faster in time and they match 
the GCS metallicity distribution at younger ages, especially 
if kch is large so metal-rich stars migrate to the Sun relatively 
rapidly. The structure of the local Hess diagram for either 
GCS stars (Fig. I16|l or Hipparcos stars (|Aumer fc BinnevI 
|2008| ) implies that the solar neighbourhood is not younger 
than 9 Gyr, so when the age is smaller than 9 Gyr the model 
is marked by a cross in Fig. 1181 Models adjacent to the 
crosses do not violate the 9 Gyr limit but are nonetheless dis- 
favoured because their local Hess diagrams yield relatively 
low likelihoods for the GCS sample. 

Four factors make it difficult to confine narrowly the re- 
quired value of kch within the part of Fig. [18] that has large 
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Figure 18. The likelihood of the GCS metallicity distribution 
in models with infall rates given by equation Upper panel: 
the values of the infall parameters /a and /b are given by the 
locations of the symbols, and the value of k^h is indicated by 
the number of sides of the polygon: 3, 4, 5, . . . for fc^h = 0.0, 
0.1, 0.2, . . .. The size of the polygon increases linearly with the 
log likelihood of the data, models with ages smaller than 9Gyr 
are marked with crosses. Lower panel: the likelihoods of models 
with /b = 0.025 and varying (/Ai'^ch)- In this panel the size of 
a symbol is a more sensitive function of likelihood than in the 
upper panel. Hexagons indicate models with best-fit ages higher 
than 12.6 Gyr, models with best-fit ages below 9.6 Gyr are crossed. 

symbols: (i) churning affects mainly the width of the metal- 
licity distribution, which has less impact on the likelihood 
than the location of its peak; (ii) churning, which is strongest 
in the inner regions of the disc, tends to saturate near the 
centre in the sense that old stars become fully shuffled; (iii) 
the GCS sample is biased against old and highly dispersed 
populations of stars, so where churning has the strongest 
effect, the observational signature is weak; (iv) the required 
churning strength is sensitive to the local metallicity gradi- 
ent which is not very well constrained by observations. 



fn our models there is no azimuthal variation in the 
metallicity of gas at a given radius, as is s uggested by recent 
observations (see lNieva fc Przvbillall200^ . which yield very 
small to negligible inhomogeneity of the fSM at a given ra- 
dius. The effect of relaxing this assumption can be gauged 
by increasing the dispersion in the measured metallicities of 
a given population of stars: if there is intrinsic dispersion 
in the metallicity of the ISM in a given annulus, the mea- 
sured metallicities of stars formed from it will reflect both 
this dispersion and measurement errors. The largest intrinsic 
dispersion in the metallicity of the ISM that would appear 
to be compatible with the data plotted in Fig.[5]is ~ 0.1 dex. 
When we combine this with measurement errors of 0.1 dex, 
we can obtain a fit to the GCS data of Fig. 1 171 that is only 
slightly worse than that provided by the standard model by 
lowering fcch from 0.25 to ~ 0.1. 

We have studied models with several values of the mass- 
loss parameter /eject and concluded that up to the largest 
values studied (/eject ~ 0.15 at i? < 3.5 kpc and 0.05 else- 
where) /eject does not have a large effect on the model's ob- 
servable properties, and is anyway degenerate with the still 
uncertain nucleosynthetic yields. However, increasing /eject 
makes it slightly easier to find an acceptable model, reflect- 
ing the fact that the yields we are using lie at the upper limit 
of the yields that are consistent with measured metallicities. 

The upper panel of Fig. [19] shows the effect on the 
fit to the GCS metallicity distribution of using a non-zero 
value of the threshold gas density, Ecrit, below which the 
SFR declines steeply. Raising Ecrit from zero to 2.5 Mq pc~^ 
changes the model prediction from the red curve of the 
standard model to the blue curve; the distribution is now 
wider and peaks at lower metallicities. The pink dotted 
curve shows the result of maximising the likelihood of the 
data subject to the constraint Ecrit ~ 2.5M0pc~'^. In this 
model /a is increased (to 0.44) and fcch is decreased (to 0.20) 
relative to the standard model. The new model provides a 
slightly worse fit to the GCS data than the standard model, 
but, as the lower panel reveals, there is a problem with us- 
ing a non-zero value of Ecrit: with minimal star-formation in 
the outer disc, the metallicity gradient of the ISM steepens 
near the edge of the star-forming regions, while further out 
the metallicity becomes constant at the intergalactic value. 
The data show no sign of this plateau, and are probably in- 
compatible with a plateau as low as [Z/'H] = — 1. Values of 
Ecrit > 2.5 Mq pc~^ are incompatible with the data because 
they bring the edge of the star-forming disc too close to the 
Sun. 

By considering the likelihoods of both the Hess diagram 
and the metallicity distribution of GCS stars, and our preju- 
dices regarding the proper value of /eject , we chose the model 
specified by Table[l]as the standard model. With this model 
the likelihood of the GCS metallicity peaks at age 11.7 Gyr. 



6 RELATION TO OTHER WORK 

Chemical evolution models of the Galaxy have a long history 
and a large literature. It would be inappropriate to attempt 
to review this literature in this section. Instead we high- 
light crucial differences with work that has most in common 
with ours, and re late our work to the analysis of the solar 
neighbourhood of iHavwoodl (|2008l V 
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Figure 19. Upper panel: the data points are the GCS counts and 
the red curve is the standard model. The broken blue curve shows 
the efEect on this model of raising S^rit from zero to 2.5 Mq pc~^. 
The pink dotted curve shows that a good fit to the data can be 
obtained for this value of Scrit- Lower panel: measurements of the 
metallicity of the ISM and the predictions of the models shown 
above. 



6.1 Comparison with earlier models 

Some of the best known models o f Gala ctic chemical evolu- 
tion are those in Chiappini et al. I (Il997l) and its successors 
IChiappini et al.l (j200ll ) and IColavitti et all (|2008h . In each 
case the disc is made up of annuli that exchange neither stars 
nor gas. Star-formatio n is driven by a K ennicutt law similar 
to equation IT}. Chia ppini et al.l (|l997t ) introduced a time- 
dependent infall rate that is superficially similar to ((6)l but 
differs from our infall rate in two important respects. First, 
the exponential with the longer time-constant is not turned 
on until 2 Gyr after the start of Galaxy formation with the 
consequence that star formation periodical ly ceases during 
the in terval 1 Gyr < t < 2 Gyr. Second, IChiappini et al.l 
(|l997t ) make the time constant &2 a linearly increasing func- 
tion of radius that va nishes at _R = 0.86 kpc (or 1.2 kpc in 
IChiappini et al.ll200ll ). whereas here 62 is constant. If we 
were to follow the prescription of Chiappini et al., our in- 
ner disc/bulge would become older, and a smaller churning 
rate would be required to bring stars more metal-rich than 
the local ISM to the solar neighbourhood. An outwards- 
increasing infall timescale enhances the metallicity gradient 
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Figure 20. Full red curve: the metallicity distribution predicted 
by the standard model for the GCS stars. Long-dashed green 
curve: the prediction of the same model for the distribution of 
the whole solar annulus. Short-dashed blue curve: the prediction 
of Chiappini et al. (1997) for the solar annulus. 



because metallicities are close to their equilibrium values, 
and these reflect the ratio of current to past SFRs. 

IChiappini et al.l l|l997h adjusted their model's free pa- 
rameters to optimise its fit to the G-dwarf metallicity distri- 
bution for stars in the s olar annulus that was determined by 
iRocha-Pinto fc Maciel 11996) from 287 stars that lie within 
25 pc of the Sun. Fig. [20] illustrates the difference between 
the metallicity distributions of stars near the Sun and in the 
entire solar annulus: the full red curve shows the standard 
model's prediction for the metallicity distribution of GCS 
stars from Fig. 1171 while the long-dashed green curve shows 
the corresponding distribution in the whole annulus. The 
distribution for the annulus is much broader than that for 
the GCS because stars with metallicities far from that of the 
lo cal ISM are typically fast-mo ving and likely to be at high 
^. iRocha-Pinto fc Macieil (|l996l ) transformed their measured 
distribution for the local sphere to the modelled globa l distri- 
bution using correction factors estimated by Sommc r-LarsenI 
(|l99ll ). which depend on the local gravitational potential 
and the velocity distributions of stars of each metallicity. 
It is clearly more satisfactory to use internally generated 
values of these distributions to predict the metallicity dis- 
tribution in the observed volume around the Sun than to 
infer the annular distribution from the measured one using 
external estimates of the velocity distributions. Moreover, 
th e short-dashed b l ue cur ve in Fig. [2S] shows the prediction 
of IChiappini et al] |l993) for the solar annulus. It declines 
much more steeply at high metallicities than the blue curve, 
predicting far too few metal-rich stars. 

The scale of the discrepancy between the blue and the 
green curves illustrat es that a model that provides an ade- 
quate fit to the data of lRocha-Pinto fc Macieil (| 19961 ) is likely 
to be incompatible with the GCS stars. This discrepancy is 
partly due to problems with the calibrat i on of the underly- 
ing dataset: iHaywoodl (|2002l ) and|Twarog (2002") pointed out 
that the metallicity calibrations of,_Schuster & Nisscn (198^) 
underlying those datasets severely underestimate the metal- 
licities of metal-rich stars. This underestimation makes the 
decline in the number of stars at [Fe/H] > steeper than it 
should be, and thus makes it easier for the data to be fitted 
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by traditional models of chemical evolution, which predict 
a sharp cutoff at high [Fe/H]. Apart from its superior cali- 
br ation, the GCS sample is 50 times larger than that used 
bv lRocha-Pinto fc Macieil (|l99Q ). so its statistical errors are 
much smaller and it is a more challenging distribution to fit. 

In the models o f IChiappini et alj ijiggTh . star formation 
ceases entirely when the surface density of gas falls below 
Ecrit, while in our models Ecrit merely marks an increase in 
the rate of decline of the SFR with decreasing gas density. 
While a cogent argument can be made for a rapid decline 
in the SFR at low gas densities , it is hard t o justify a dis- 
continuity in the rate. IChiappin i et al.l l|200lh conclude that 
a non-negligible value of Ecrit plays an essential role in fit- 
ting the data. In fact, they set Ecrit = 7M0pc~^, with the 
consequence that star formation in the solar neighbourhood 
was constantly stopping and starting, both during the first 
2 Gyr and the last 4 Gyr of Galactic history (see their Fig. 4). 
Since the metallicity of the ISM declines when star forma- 
tion has been switched off, this erratic behaviour broadens 
the stellar metallicity distribution. In our models the SFR is 
steady and in fact setting Ecrit ~ 7Mq pc~^ would result in 
almost no stars forming at i? > 10 kpc because the flow of gas 
through the disc would prevent the surface density building 
up to 7M0pc~^. Hence in our models star-formation can 
occur at large radii only if Ecrit is set to a small or van- 
ishing value. An important d i fferen ce between our models 
and those of iGhiappini etall l|200ll ) is that in our models 
the SFR is a smooth function of time and the bimodality in 
[g/ Fe] is achieved w it hout a dip in the star-formation rate. 

IChiappini et al.l |l993) give for the inner and outer 
parts of their models the radial gradients in the abundances 
of several elements. These gradients are largest in the inner 
regions, but even there they are much smaller than in our 
models: for example at 12 Gyr the inner gradients of [0/H] 
and [Fe/H] are —0.023 and — 0.027 dexkpc"^ compared to 
values ~ — 0.08dexkpc~^ and ~ — 0.11 dexkpc"^ obtained 
here. Our larger gradients are a direct consequence of the 
advection inwards of the products of nucl eosynthesis. The 
GCS s tars show a similar gradient in [Z/H] : lHolmberg et al.l 
l|2007l) derive a gradient o f — 0.09dexkpc~^. 

iGolavitti et all (|2008l ) used infall rates measured from 
simulations of clustering cold dark mat t er (C DM) in 
slightly updated models of IChiappini et al] (|l997l ). No em- 
pirically determined infall rate gave such satisfactory re- 
sults as the earlier double-exponential rate. The models 
were again compared to th e metallicity distribution of 
iRocha-Pinto fc Macieil (|l996l ) and in most cases provided 
inadequate fits to the data. The empirically measured in- 
fall rates are very irregular in time, with the result that the 
model experiences powerful bursts of star formation. These 
lead to large excursions in the predicted plots of [O/Fe] ver- 
sus [Fe/H] for which there is no evidence in the data. More- 
over, studies of the past SFR in the disc show no signs of ma- 
jor bursts of star formation a few gigayears ago. The disap- 
pointing results of this study suggest that the rate at which 
gas joins the disc is not simply the ratio of masses of baryons 
and dark-matter times the dark-matter accretion rate. In 
fact, the wealth of evidence that t he majority of baryons 
are still in the inter galactic medium (|Persic fc Salucci|[l993 : 
iFukigita et al.|[l998l ) is a clear indication that galaxies do not 
acquire gas as fast as this naive calculation suggests; rather 
gas is stored in the warm-hot intergalactic medium (WHIM) 



and from there accreted at a still uncertain rate. It seems 
unlikely that chemical evolution models will be successfully 
coupled to cosmological clustering simulations until we have 
understood the complex interface between the WHIM, cold 
infa ll and galactic fountains . 

iNaab fc Ostrikerl (|2006l ) determined the infall rate by 
assuming that the surface density of the disc is always ex- 
ponential, but with a scale length that is proportional to 
Vc/H , where Vc is the Galaxy's circular speed (taken to be 
constant) and H{t) is the Hubble parameter. The infall rate 
at all times and places was fixed by assuming that the cen- 
tral surface density of the disc is constant at its current 
value. The model's observables were calculated by using the 
Kennicutt law ([T} to convert gas to stars. A very simple 
approach to chemical evolution was employed, in which ma- 
terial does not move between annuli and only the overall 
metal content Z was followed. In these models the metallic- 
ity gradient in the ISM becomes less steep over time, and is 
now ~ — 0.04dexkpc~^. Their model predicts for the solar 
neighbourhood fewer metal-rich stars and more metal-poor 
stars than are in the GCS. 

In all these models, solar-neighbourhood stars should 
satisfy a well defined metallicity-age relation, and the G- 
dwarf metallicity distribution is simply the result of com- 
bining this relation with the SFR-age relation. As we have 
seen, the observed width of the local metallicity distribu- 
tion is approximated by exploiting irregular time evolution 
of th e metallicity of the local ISM (cf. also IPortinari et al.l 
Il998h . As Fig. [5] illustrates, our models solve this problem 
in an entirely different way, and using a simpler history of 
star formation. 

6.2 Haywood's analysis of the solar 
neighbourhood 

iHavwoodl l|2006l 'l and especially iHavwoodl l|2008l ') has crit- 
ically re-examined the age-metallicity distribution of the 
GCS stars and concluded that the data are only consistent 
with the existence of a well defined age-metallicity relation 
for disc stars younger than ~ 3 Gyr; for such young stars the 
spread in [Fe/H] is consistent with the expected dispersion 
in the metallicity of interstellar gas at a given radius. At ages 
larger than 3 Gyr the width of the metallicity distribution 
is larger than can be accounted for by measurement errors 
and inhomogeneity of the ISM. In particular, there are stars 
with ages ~ 5 Gyr that have [Fe/H] — 0.5, and stars with 
ages < 7 Gyr that have [Fe/H] = —0.5. The older the age 
bracket that one examines, the wider the range of metal- 
licities present. Our predicted age-metallicity distribution 
(Fig. 16} is in good agreement with that derived bv IHavwoodl 
(200^ 

iHavwoodl l|2006l ) finds that when his revised ages are 
used for GCS stars, the metallicities of thick disc stars in- 
crease as their ages decrease. That is, these stars point 
to rapid self-enrichment of the thick disc. In fact, Hay- 
wood argues that the thick disc is not the relic of some 
captured satellite(s) but an integral part of the Galaxy's 
disc and has played a central role in the chemical evolu- 
tion of the thin disc. Chemical evolution models should 
treat the disc as a whole, not just individual parts. Our 
results strongly underline this conclusion from a theoretical 
perspective. The metallicity-age plot shown in Fig. 1(b) of 
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iHavwoodl (|2008l ) is satisfyingly similar to the lower panel 
of our Fig. |6l in both figures the stellar density is highest 
around (r — 2 Gyr, [Fe/H] = —0.1) and in this region the 
ridge line gradually drops to the right. At older ages the dis- 
tributions become broader, being confined by [Fe/If] ~ 0.4 
and —0.5. At ages greater than 10 Gyr both distributions 
reach down to [Fe/H] — —1. In fact the small differences be- 
tween our figure and that of Haywood are readily accounted 
for by the substantial errors in measured stellar ages. This 
fit is remarkable because the model parameters were chosen 
without reference to measured stellar ages. 

An argument sometimes advanced for a dichotomy be- 
tween the thick and thin discs is the existence of two se- 
quences in the ([a/Fe], [Fe/H]) plane (Fig. H]). This dia- 
gram suggests that the last thick-disc stars to form had 
higher abundances than the first thin-disc stars to form. 
To explain this finding in the context of a conventional 
chemical-evolution model, sudden dilu tion of the ISM b ; 
a massive gas-rich accretion is required (iBensbv et al.|[200 
iReddv et aLlbOod ). Against this proposal HavwoodI (|200i 
objects that the bulk of the oldest thin-disc stars have 
[Fe/H] ~ —0.2 and there is no evidence that the most 
metal-poor thin-disc stars are particularly old. Our mod- 
els show that the observed structure of the ([a/Fe], [Fe/H]) 
plane ari s es na turally when radial migration is allowed. 
IHavwoodl l|2008l l examined the orbital parameters of stars 
of various metallicities and showed that local thin-disc stars 
with metallicities that overlap the metallicity range of the 
thick disc have higher angular momenta than more typical 
thin-disc stars. Similarly, he found that stars in the high- 
metallicity tail of the local metallicity distribution have low 
angular momenta. Even though churning could in principle 
eradicate the correlation between angular momentum and 
metallicity, our models reproduce these correlations (Figs [9] 
and llOp because blurring makes a sufficiently large contribu- 
tion to bringing these chemically anomalous stars near the 
Sun.^ 

llvezic et al.l (|2008l ) also argue that in the SDSS data the 
kinematic properties of the thick disc evolve continuously 
with distance from the plane in a way that suggests that 
the t hick disc joi n s con tinuously to the thin disc. By con- 
trast, IVeltz et ai] (|2008l ') argued for a clean break between 
the thin and thick discs on the basis of shallow local minima 
in the density of 2MASS stars n{z) in seen at t he Galactic 
poles as a function of photometric distance z. As lVeltz et al.l 
l|2008ll show, minima in n{z) are not expected if the disc is 
a superposition of two exponential structures, but the min- 
ima yield a clean discontinuity in the distribution of velocity 
dispersions when multi-component isothermal distribution 
functions are used to model the data. The very unexpect- 
edness of the minima makes the modelled break extremely 
clean. It will be interesting to see whether the distribution 
of measured radial velocities of stars in the RAVE survey 
substantiate these model velocity distributions. 



7 CONCLUSIONS 

It is now more than forty years since the theory of stellar 
evolution attained the level at which it became possible to 
model the chemical enrichment of the ISM. From the be- 
ginning of that endeavour measurements of the abundances 



of individual solar-neighbourhood stars have played a key 
role because a star preserves like a time capsule the state of 
the ISM at the remote epoch of its formation. Considerable 
theoretical and observational efforts have been devoted to 
probing the history of t he Galaxy with this connection. 

Half a century ago, iRoma nl (l950l . [l953 ) and others dis- 
covered the connection between the kinematics and chem- 
istry of stars, yet curiously little has been done to include 
kinematics in models of chemical evolution. The general pre- 
sumption has been that each annulus of the disc evolves in- 
dependently of others, and the well-known correlations be- 
tween chemistry and kinematics can be understood as aris- 
ing through the stochastic acceleration of stars: older stars 
tend to have larger random velocities and lower metallicities. 
No effort was made to develop greater diagnostic power by 
simultaneously modelling chemistry and kinematics. 

The continued use of mutually independent annuli by 
modellers of chemical evolution is surprising given that it 
was from the outset recognised that many stars are on signif- 
icantly non-circular orbits that each radial half-period cover 
more than a kiloparsec in radius. In fact, it has generally 
been assumed that the radial velocity dispersion within the 
disc rises as one moves inward to values that lead to ra- 
dial excursions of several kiloparsecs (see Fig. |4|. Moreover, 
observations have long indicated th at galactic discs have sig- 
nificant me tallicity gradients (e.g. IVila-Costas fc Edmunda 
Il992l : lde JonE.1996 ). so radial migration of stars is bound 
to leave a signature on the metallicity distribution of solar- 
neighbourhood stars. We have called this aspect of radial 
migration "blurring" . 

The present study ow es its impetus to the discovery by 
ISellwood fc BinnevI (|2002D that radial migration is a more 
potent process than mere blurring: the dominant effect of 
transient spiral arms is not to heat the stellar disc as had 
been supposed, but to cause stars either side of corotation 
to cha nge places without moving to eccentric orbits ( "churn- 
ing" ) . ISellwood fc BinnevI l|2002l) did not demonstrate that 
gas participates in churnin g, but they argued th at it must 
on dynamical grounds, and iRoskar et af] (20083) found ev- 
idence that this was the case in their N-body-SPH simu- 
lations of galaxy formation. Since churning is an aspect of 
spiral structure that can only be probed through its impact 
on chemical evolution, we wanted chemical-evolution mod- 
els that included churning, and logically it was natural to 
extend these models to include both blurring and radial gas 
flows. 

An ingredient of our models that might be controversial 
is the introduction of radial gas flows. We have absolutely 
no reason to expect that the infall proflle is exponential, so 
if discs are to be exponential this must be the result of flows 
through the discs redistributing mass within the disc. In fact, 
as gas streams through spiral arms it dissipates energy in 
shocks that is ultimately gravitational energy that becomes 
free as the gas surrenders angular momentum to the stars 
and drifts in wards. Hence at so me level inward gas flows are 
mandatory l|Lacev fc Fal]||l985l ). 

Unfortunately, the theory of galaxy formation has yet 
to advance to the point at which it can prescribe the spatial 
and temporal structure of gas accretion, so it is necessary 
to parametrise accretion in some way. The accretion process 
must be constrained to result in the formation of the ob- 
served stellar and gaseous discs. Our accretion Scheme AB 
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satisfies this constraint for all values of the parameters, but 
it is inevitably acausal in that the formation mechanism is 
being driven by its known outcome. While its acausality is 
unattractive, Scheme AB is a fiexible parametrisation that 
enables us to form exponential discs for a variety of differ- 
ent assumptions about the radial density of infalling gas, and 
the resulting radial profiles of infall and gas-fiow (Fig. [2]) are 
entirely plausible. 

The impact that radial migration has on the local metal- 
licity distribution obviously varies with the magnitude of the 
metallicity gradient in the ISM, which in turn depends on 
the gas fiow within the disc and therefore the radial infall 
profile. For this reason the most important parameters of 
our models are /a and /b, which control the distribution of 
infalling gas. 

The models provide good fits to the GCS counts of stars 
as functions of [Fe/H], Mv, Tcff and stellar age, as well as 
reproducing the correlations between tangential ve locity and 
abund ance patterns that have been pointed out bv lHavwoodI 
l|2008h . These fits are achieved by churning, blurring and 
radial fiows working together. They depend on the existence 
of an appreciable metallicity gradient in the ISM, which is 
established by the radial flow of gas, and they depend on 
radial mixing of stars by blurring and churning. The steeper 
the metallicity gradient in the gas, the smaller the effect 
of churning can be, but for any observationally consistent 
metallicity gradient, churning has a non-negligible role to 
play. 

The models describe the coevolution of the thick and 
thin discs, and presume that thick-disc stars were formed 
in the Galaxy rather than accreted from outside. Given 
the simplicity of our assumptions, the extent to which a 
dichotomy between an a-enhanced thick disc and a solar- 
type thin disc automatically manifests itself in the mod- 
els is remarkable. In particular, in the solar annulus the 
distributions of [O/H] at given [Fe/H] are bimodal in the 
range of [Fe/H] associated with the overlap of the two discs 
(Fig. the vertical density profile can be represented at 
the sum of two exponentials, and at _R < _Ro the radial 
density profile becomes fiatter at lager distances from the 
plane (Fig.jH]). None of these characteristics is dependent on 
our choice of a double exponential for the time dependence 
of infall: a model in which the gaseous mass, and therefore 
star-formation rate, is held constant also has these features. 
They are consequences of the ~ 1 Gyr timescale of type la 
supernovae and the secular heating and churning of the disc. 

The models assume that scattering of stars increases 
the velocity dispersion of a coeval population as t^^^ , but 
the models go on to predict that within the solar neighbour- 
hood velocity dispersion increases as a higher power of age, 
roughly f"'^'' in agreement with what is found from Hippar- 
cos stars with good parallaxes. This finding may reconcile 
scattering theory, which cannot readily explain an exponent 
in excess of 1/3, with observations. The key point is that 
radial mixing brings to the solar neighbourhood stars born 
at small radii, where the velocity dispersion is undoubtedly 
large. 

The nucleosynthetic yields from each generation of stars 
are still significantly uncertain. Our philosophy has been to 
use standard values from the literature rather than exploit 
uncertainties in the yields to tune the models to the data. 
The yields we are using are in the upper region of those that 



can provide adequate fits to the data, with the consequence 
that increasing the value of the mass-loss parameter /eject, 
which is degenerate with the magnitude of the yields, makes 
it easier to fit the data. Our yields are surely not exactly 
right and consequently some of the properties the models 
derive from them will be in error. On account of these un- 
certainties we have suppressed the predictions of our code for 
the abundances of certain elements, most notably carbon. 

The discovery that three-dimensional, non-equilibrium 
models of the solar atmosphere require the metal abundance 
of the Sun to be Zq = 0.012 - 0.014 jOrevesse et al.ll2007h 
poses a major problem for this field. A prerequisite for suc- 
cessful chemical modelling is a consistent metallicity scale 
for both stars and the ISM. At present the only consistent 
scale is the traditional one on which Zq — 0.019, so this is 
the one we have used. 

If a new scale were established on which all metallicities 
were significantly lower, viable models could be produced 
by lowering the yields. A straightforward way to do this 
would be to lower the maximum mass in the IMF: reducing 
this mass from 100 M© to 50 Mq would reduce yields by 
~ 30 percent in line with the proposed reduction in Zq. 
The oxygen yield would come down fastest, reducing [0/Fe] 
by ~ 0.2 dex. 

Although we believe that this study represents a signif- 
icant advance on all previous models of Galactic chemical 
evolution, it is highly imperfect. Some major weaknesses of 
our work are the following. 

1. We have assumed that the probability of mass inter- 
change between rings is proportional to the product of the 
rings' masses. This probability reflects the number and in- 
tensity of spiral feature with corotation at that radius, and 
should be a function of both mass and velocity dispersion. 
A further study of self-gr avitating discs similar to that of 
ISellwood fc BinnevI (|2002D would be necessary to determine 
this function. 

2. We have assumed that the vertical and radial motions 
of stars decouple. This assumption has a significant impact 
on both the relation between age and velocity dispersion 
in the solar neighbourhood and on the predicted vertical 
density profile in the solar annulus, which is interesting in 
itself and impacts on the selection function of GCS stars 
and thus on our choice of standard model. The assumption 
is unjustifiable for stars on eccentric orbits, which do play an 
important role in the model fits. Unfortunately, a sounder 
treatment is impossible until a better approximation to the 
third integral of Galactic dynamics is available. It is our 
inte ntion to resolve this proble m through "torus modelling" 
(e.g. iMcMillan fc BinnevI [ioosl ). 

3. Our models include radial mixing of gas through 
churning and viscous inspiralling, but do not include ra- 
dial redistribution of g as by the Galactic fountain (e.g. 
iBeni amin &: Danlvl[l997l ) . A significant body of evidence in- 
dicates that star formation drives neutral hydrogen to kilo- 
parsec heights above the plane. NGC891, which is not dis- 
similar to the Galaxy, has ~ 25 per c ent of its neutral hydro - 
gen more than 1 kpc from the plane (jOosterloo et al.ll2007l ). 
This extraplanar gas must move over the plane on nearly bal- 
listic trajectories, and in the absence of interaction with the 
corona (gas at the Galaxy's virial temperature ~ 2 x 10^ K), 
the gas must return to the plane further out than its point 
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of ejection. However, observations suggest that neutral gas 
above the plane is actually flowing inward rather than out- 
ward, presumable as a result of interaction with the corona 
l|Fraternali fc Binnevlliooa ') . The mass of extraplanar gas is 
so large and the timescale for its return to the plane so short 
that whichever way this gas flows, it has a considerable po- 
tential for radially redistributing metals. Extraplanar gas 
must be ejected from the disc by supernova-heated gas that 
is probably highly metal-rich. Some of this gas will be lost 
to the Galaxy as we have assumed, but some of it will return 
to the plane with infalling gas. Again there is potential here 
for significant radial redistribution of metals that has been 
ignored in the present study. 

4. Our treatment of the inner Galaxy is unacceptably 
crude in that we have replaced the bar/bulge with a disc. 
Unfortunately introducing the bar opens a Pandora's box of 
complexities, and at the present time is probably only feasi- 
ble in the context of a part icle-b ased model such as th ose of 
ISamland fc GerhardI (|2003l ) and lRoskar et al.l (|2008bD . Our 
hope is that our imaginary central disc has a similar impact 
on the chemodynamical state of the local disc to the com- 
bined impact of the giant star- forming ring at i? ~ 4kpc, the 
gas-deficient region interior to this ring and the star-forming 
X2 disc at i?<200pc. 

From this list of shortcomings of our models it is clear 
that we are still far from a definitive account of the Galaxy's 
chemical evolution. We shall be satisfied if we have con- 
vinced the reader that the interplay between dynamics and 
chemistry is so tight at to be indissoluble. This fact is at one 
level inconvenient because it undermines the value of conclu- 
sions drawn from traditional models of chemical evolution. 
But at another level it represents an opportunity to learn 
more. The connection between the kinematics of stars and 
the compositions of stars and gas, which can be measured 
in great detail, involves three areas about which we are too 
ignorant: the distribution of dark matter within and around 
the Galaxy, the Galaxy's history of assembly, and the nature 
of the Local Group's IGM. By building dynamical models 
of the Galaxy that have chemical evolution built in to the 
basic structure, we should be able to make decisive progress 
with one of the major problems of contemporary astronomy. 
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Figure 21. The distribution of a abundances predicted by equa- 
tion |[2T||. 



a coeval group of stars formed at t' generates a rate of type la 
supernovae that vanishes for t < t' + to and is subsequently 
oc e"''"''^' Then given that the rate of core-collapse SNe is 
proportional to the SFR, the rate of production of Fe is 
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where 6 and c are constants. Integrating this production 
rate, we obtain the iron mass at time t > to as 
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In the approximation that SNIa do not contribute a ele- 
ments, and that the delay in the production of these ele- 
ments by a stellar population is ^ 1/fc, the mass of a ele- 
ments is Ma = (1 — e~^^)a/K where a is a constant. Then 
equation (|21|l predicts that a = Mc/Myc equals a/h for 
t < to and then drops rapidly towards its asymptotic value, 

.(oo) = A + ^e-o. (22) 

Fig. [21] plots the distribution at t = 13 Gyr of stars over 
Q when the initial and asymptotic values of a are set to 
0.4 and 0.2 and the other parameters are K — 1/7 Gyr, 
k — 1/1.5 Gyr and to = 0.3 Gyr, which allows 0.15 Gyr for 
white dwarfs to form and 0.15 Gyr for them to accrete prior 
to deflagrating. 



APPENDIX: ORIGIN OF BIMODAL [O/FE] 
DISTRIBUTIONS 

We provide an analytic model of the development of the 
bimodal distributions of [O/Fe] evident in Fig. [5] We assume 
that star formation starts at t = and that the SFR is oc 
e~^* for t > 0. Consistent with equation ([5]), we assume that 



